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Abstract 



> . 

, The effect of shadowing on the early state of ultrarelativistic heavy ion col- 

lisions is investigated along with transverse energy and hard process produc- 
■ tion, specifically Drell-Yan, J/tp, and T production. We choose several parton 

distributions and parameterizations of nuclear shadowing, as well as the spa- 



tial dependence of shadowing, to study the influence of shadowing on relevant 
observables. Results are presented for Au+Au collisions at ^/snn = 200 GeV 



Q_i| and Pb+Pb collisions at ^snn = 5.5 TeV. 

43 . I. INTRODUCTION 

> '■ 

Experiments have shown that the proton and neutron structure functions are modified 
in the nuclear environment. The modification depends on the parton momentum fraction x. 
For medium x, 0.3 < x < 0.7, the nuclear parton distributions are depleted relative to those 
in isolated nucleons. For intermediate x, 0.1 < x < 0.3, the distributions are enhanced, an 
effect known as antishadowing. Finally, for small x, x < 0.1, the nuclear depletion returns. 
We refer to the entire characteristic modification as a function of x as shadowing. To date, 
most measurements of shadowing have studied charged partons, quarks and antiquarks, 
through deep-inelastic scattering (DIS), eA — > e'X, while the behavior of the nuclear gluon 
distribution has been inferred from the modifications to the charged partons. 

Almost all of these measurements were blind to the position of the nucleons in the nucleus. 
However, most models of shadowing predict that the structure function modifications should 
be correlated with the local nuclear density. For example, if shadowing is due to gluon 
recombination, it should be proportional to the local nuclear density. The only experimental 
study of spatial dependence of parton distributions relied on dark tracks in emulsion to tag 
more central collisions M. They found evidence of a spatial dependence but could not 
determine the form. 
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This paper studies the effect of shadowing and its position dependence in ultrarelativistic 
Au+Au collisions at a center of mass energy of 200 GeV per nucleon, as will be studied at 
the Brookhaven Relativistic Heavy Ion Collider (RHIC) and in 5.5 TeV per nucleon Pb+Pb 
collisions expected at the CERN Large Hadron Collider (LHC). We determine the initial 
quark, antiquark and gluon production and the first two Et moments of minijet production 
for two commonly used parton distributions, with three shadowing parameterizations. Fol- 
lowing previous calculations, we find the initial energy density and the average energy per 
particle. We critically examine the concept of fast thermalization in these collisions. 

The spatial dependence of shadowing is reflected in particle production as a function of 
impact parameter, b, which may be inferred from the total transverse energy, Et, produced 
in a heavy ion collision We discuss the relationship between Et and b including both 
hard and soft contributions. We then consider the effect of shadowing on the production 
of hard probes such as J/ip, T, and Drell-Yan dileptons as a function of b. These latter 
calculations complement studies of shadowing in open charm and bottom production || . 

Section 2 discusses the initial state nuclear parton distributions, including shadowing 
and its spatial dependence. Section 3 then considers minijet production and the effect on 
initial conditions for further evolution of the system. Section 4 is devoted to the relationship 
between transverse energy and impact parameter. Section 5 discusses J/ip and Drell-Yan 
production and their sensitivity to the nuclear parton distributions. Finally, section 6 gives 
some conclusions. 



II. NUCLEAR PARTON DISTRIBUTIONS 

The nuclear parton densities, F Z A , are assumed be the product of the nucleon den- 
sity in the nucleus Pa( s ), the nucleon parton density f^(x,Q 2 ), and a shadowing function 
S l (A, x, Q 2 , r, z) where A is the atomic mass number, x is the parton momentum fraction, 
Q 2 is the interaction scale, and r and z are the transverse and longitudinal location of the 
parton in position space with s = \J\r\ 2 + z 2 so that 

F A (x, Q 2 , r, z) = p A (s)S\A, x, Q 2 , r, z)f l N (x, Q 2 ) . (1) 

In the absence of nuclear modifications, S i (A, x, Q 2 , r, z) = 1. The density of nucleons in 
the nucleus is given by the Woods-Saxon distribution 

l + u(s/R A ) 2 

Pa{s) = Po— r, w7/T ' 2 

1 + exp[(s - RA)/d\ 

where the nuclear radius, Ra, skin thickness, d, and oblateness, to, are determined from 
low energy electron-nuclear scattering ||. The central density is determined by the nor- 
malization / d 2 rdzpA{s) = A. Results are given for Au+Au collisions at RHIC and Pb+Pb 
collisions at the LHC with R^ u = 6.38 fm and Rp^ = 6.62 fm respectively. 

The densities of parton i in the nucleon are obtained from fits to DIS data. These fits 
are necessary because the distributions at the initial scale Qo are nonperturbative. However, 
the parameterizations of f± are only reliable where measurements exist. The continually 
improving DIS data from HERA shows that uncertainties still exist at small x. Therefore, 
we consider two different parton distribution sets. Both are chosen because they are leading 
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order, LO, sets which are more consistent with a leading order calculation. They also have 
a relatively low initial scale. The GRV 94 LO || distributions have a lower scale, Qq = 0.63 
GeV, than the MRST LO |J central gluon distribution with Q = 1 GeV. Figure [I] compares 
the valence, qy = uy + d v , sea quark, q s = 2(u + d + s), and the gluon, g distributions 
at Q 2 = 4 GeV 2 of the two sets. In the low x region probed by the LHC the valence and 
sea quark distributions in both sets are similar. However, the MRST LO gluon distribution 
is less than half as large as the GRV 94 LO gluon distribution. As we show in the next 
section, the gluons dominate particle production at the LHC. Thus the low x gluon density 
will significantly affect the initial conditions obtained for these high energies. At x ~ 0.01 
the parton densities are well known so that the two sets are similar and the choice of the 
proton parton distributions do not strongly influence the initial conditions at RHIC. 

Shadowing is an area of intense study with numerous models available in the literature 
|J • However, none of the models can satisfactorily explain the behavior of the nuclear parton 
distributions over the entire x and Q 2 range. Therefore, we choose to use parameterizations 
of shadowing based on data. We use three different fits, all based on nuclear DIS data. As 
in DIS with protons, the nuclear gluon distribution is not directly measured and can only 
be inferred from conservation rules. The first parameterization, Si(A, x), treats quarks, 
antiquarks, and gluons identically without Q 2 evolution |TIJ. The other two evolve with 



Q and conserve baryon number and total momentum. The S\{A, x, Q ) parameterization, 
starting from the Duke-Owens parton densities modifies the valence quarks, sea quarks 



and gluons separately and includes Q 2 evolution for Q = 2 < Q < 10 GeV Hl2| . The third 
parameterization, S\(A, x, Q 2 ), is based on the GRV LO |13| parton densities. Each parton 



type is evolved separately above Qo = 1.5 GeV [n,D|]. The initial gluon distribution in 
63 shows significant antishadowing for 0.1 < x < 0.3 while the sea quark distributions are 
shadowed. In contrast, S 2 has less gluon antishadowing and essentially no sea quark effect 
in the same x region. Since S3 includes the most recent nuclear DIS data, it should perhaps 
be favored. Figure compares Si, 5*2 and S3 for Q = Q Q and Q = 10 GeV. 

The remaining ingredient is the spatial dependence of the shadowing. Unfortunately, 
there is little relevant data. Fermilab experiment E745 studied the spatial distribution of 
nuclear structure functions with vN interactions in emulsion. The presence of one or more 
dark tracks from slow protons is used to infer a more central interaction [[| . For events with 
no dark tracks, no shadowing is observed while for events with dark tracks, shadowing is 
enhanced over spatially independent measurements from other experiments. Unfortunately, 
this data is too limited to be used in a fit of the spatial dependence. 

Most models of shadowing predict that the nuclear parton densities should depend on 
the interaction point within the nucleus. In one model, at high parton density gluons and 
sea quarks from one nucleon can interact with partons in an adjacent nucleon |16| so that 
shadowing is proportional to the local density, Eq. @ |||5|]. Then 

S l ws = S l (A,x,Q 2 ,r,z) = l + N ws [S l (A,x,Q 2 ) - 1]^ , (3) 

Po 

where A" ws is chosen so that (A/A) J d 2 rdzpA(s)S^ s = S\ At large distances, s 3> Ra, the 
nucleons behave as free particles while in the center of the nucleus, the modifications are 
larger than the average value S\ 

In another approach, shadowing stems from multiple interactions by an incident parton 



Tq . Parton-parton interactions are spread longitudinally over a distance known as the 
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coherence length, l c = l/2mx, where m is the nucleon mass For x < 0.016, l c is greater 
than any nuclear radius and the interaction of the incoming parton is delocalized over the 
entire trajectory. The incident parton interacts coherently with all of the target partons 
along this interaction length. At large x, l c <C R A and shadowing is proportional to the 
local density at the interaction point, while for small x, it depends on the density integrated 
over the incident parton trajectory. Both formulations reproduce the spatially independent 
shadowing data quite well. Unfortunately, the available data |2[] is inadequate to test these 
theories. 

Because of the difficulty of matching the shadowing at large and small x while maintaining 
baryon number and momentum conservation, we do not include the multiple scattering 
model explicitly in our calculations. However, we do consider the small x and large x limits 
separately. Equation (|3p corresponds to the large x limit. In the small x regime, the spatial 
dependence may be parameterized 

S p {A,x,Q 2 ,r,z) = l + N p (S l (A,x,Q 2 )-l) Jdzp A (f,z). (4) 

The integral over z includes the material traversed by the incident nucleon. The normaliza- 
tion requires (I /A) J d 2 rdzp A (s)S' p = S\ We find N p > A ws . 

There are a number of difficulties with the coherent-interaction picture. While traversing 
the formation length, both the incident and the produced partons will undergo multiple 
interactions, which will reduce the effective coherence length, analogous to the Landau- 
Pomeranchuk-Migdal effect [|19|]. Also, the picture of a single incident parton interacting 
with a static nucleus is inappropriate in heavy ion collisions since the parton density rises 
rapidly as many interactions occur simultaneously. A step-by-step calculation cannot solve 
this problem because non-local depictions of heavy ion collisions are inevitably Lorentz 



frame dependent [20]. Finally, in a model where the parton densities are spread out over an 
x-dependent distance, baryon number is not locally conserved. 

We previously considered a variant, S R , where shadowing is proportional to the thickness 
of a spherical nucleus at the collision point ||, 



itlA^.M = 1+ *^»,<n - W 1 - . (5) 

I 1 r>R A 

The normalization, Ar, obtained after averaging over p A (s), is similar to N p . This model 
suffers from a discontinuous derivative at r = R A with no shadowing predicted for r > R A , 
but is otherwise fairly similar to S p . 

Figure ^| compares the radial dependence of 5ws> S p , and Sr for S l (A, x, Q 2 ) = 0.7. 
For the comparison, S\ys is evaluated at z = 0. The S p and Sr results are very similar 
except near the nuclear surface where they differ by ~ 10%. Later we compare calculations 
of the first moment with S\ys an d S p and show that the two results are very similar. 
Calculations using S p and Sr would be in closer agreement, effectively indistinguishable. 

Other mechanisms such as nuclear binding have also been suggested as possible explana- 



tions of shadowing |ZTJ. These calculations can explain only a small fraction of the observed 



effect p2 |, at least for x > 0.1. However, many of these models would also predict some 



spatial dependence. 
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Given the difficulties of matching spatial dependencies for different x and A while pre- 
serving baryon and momentum conservation in the multiple interaction model, we focus our 
calculations on the local density model, and perform most of our calculations using Sws- 
However, as we will show, the calculations are relatively insensitive to the exact parame- 
terization suggesting that heavy ion collision studies will not distinguish between different 
models. 

For simplicity, we will refer to homogeneous (without spatial dependence) and inhomo- 
geneous (position dependent) shadowing. 



III. INITIAL CONDITIONS IN A + A COLLISIONS 



At RHIC and LHC perturbative QCD processes are expected to be an important compo- 
nent of the total particle production. At early times, n ~ 1/pr < 1/Po ~ 0.1 fm for po ~ 2 
GeV, semihard production of low pt minijets will set the stage for further evolution 
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Copious minijet production, especially gluonic minijets, in the initial NN collisions has been 
suggested as a mechanism for rapid thermalization, particularly at the LHC. We critically 
examine this idea, with special attention to the effects of shadowing on these expectations. 

Minijet production is calculated from the jet cross section for px > po- At leading order 
the rapidity, y, distribution of a parton flavor / produced in the parton subprocess ij — > kl 
in AB collisions is 12T 



d%di°dy = K]et I d PTdy2dzdz'^2x 1 F l A (x 1 ,PT,r,z)x 2 Ff(x 2 ,PT,b- f,z') 

(ii) 



x 



1 + 5, 



(6) 



-P 2 t( 



where t = —p\{l + e ^ y2 - ) ) and u 
and y 2 are p\ < p\ < sattv/(4 cosh 2 y) and m(r. 
ln(r Po + 



1 + e 



{y-yzY 



-v 



PT 



y/s NN /p T andr po = y/s NN /2p Q 



. The limits of integration on p\ 
) < V2 < l n ( r PT — e v ) where \y\ < 
The sum over initial states includes 



all combinations of two parton species with three flavors while the final state includes all pairs 
without a mutual exchange and four flavors (including charm) so that a s {pT) is calculated 
at one loop with four flavors. The factor 1/(1 + 5m) accounts for the identical particles in 
the final state. The factor K ]et in Eq. is the ratio of the NLO to LO jet cross sections 



Previous analysis of high p T jets predicted 
A more recent NLO calculation of minijet production 
found K- Jct pa 2 at both RHIC and LHC ||26|| . Assuming K )et = 1, as in Ref. p4| , gives a 



and indicates the size of the NLO corrections 
1.5 at LHC energies 



jet 



ID- 



conservative lower limit on minijet production. The cutoff po, represents the lowest pt scale 
at which perturbative QCD is valid. There is some uncertainty in the exact value of po 

However, 2 GeV should be a safe value for 
. The effects of different choices for p will 



which can be constrained by soft physics |27| 



heavy ion collisions, especially at the LHC |24 
be discussed later. 

The parton densities are evaluated at scale pt, with x values at as low as x\ >2 ~ 
2po/\/snn ~ 7 x 1CT 4 at y = y 2 = in Pb+Pb collisions at 5.5 TeV/nucleon. At higher 
rapidities, x± or x 2 can be even smaller. Thus the small x behavior of the parton densities 
strongly influences the initial conditions of the minijet system. 
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The resulting minijet rapidity distributions are shown in Figs. ^-0 for the two sets of 
parton distributions at the LHC and RHIC both without shadowing and with homogeneous 
shadowing. Shadowing can reduce the number of produced partons by up to a factor of two 
at the LHC, depending on the parameterization and the parton type. The smallest effect is 
observed with the newer S3 parameterization. At the lower RHIC energy, xxp ~ 0.02, and 
shadowing is smaller, as is shown in Fig. |2|. Due to the strong antishadowing, gluons are 
actually enhanced with 5*3. 

Since each collision has two final state partons, the total number of partons of flavor / 
at impact parameter b is 

^(6,p„) = 2^ (7) 

where da? (po) j d 2 b is the integral of Eq. over d 2 r and dy normalized so that 

d<J f (p ) n da f (p ^ 



d*bdy dV = 2 ^b~ (8) 

because there are two final-state partons in each collision. The total hard scattering cross 
section as a function of impact parameter is the sum over all parton flavors so that 

When S — 1 or S = S' l (A, x, Q 2 ) the spatial dependence factorizes, the per nucleon cross sec- 
tion is independent of b, and the total cross section scales with the nuclear overlap function, 
TAB^b), p8| . The overlap function is the convolution of the nuclear density distributions f| 



T AB (b) = J d 2 fT A (f)T B (b-f) (10) 

with the nuclear thickness function Ta{t) = J dzpA^z^r). For AA collisions, 7^4(0) ~ 
A 2 /(tiR a ) oc A 4 / 3 . The transverse area of the system and the initial volume at b = are 

A T = tiR 2 a (11) 
Vi = A T Ayr, = A T Ay/p , (12) 

where Tj = 1/po and Ay is the rapidity range. 

Parton production saturates when the transverse area occupied by the partons is larger 
than the total transverse area available. The total number of partons produced in the 
collision is the sum over flavors, 

A^H^iV^po) (13) 
/ 

In a b = collision, the partons occupy a transverse area ttN (0) jp\. Saturation occurs when 
the area occupied by partons is equivalent to the transverse area of the target in a symmetric 
heavy ion collision at b = 0, iV (0) > R\Po- I n Pb + Pb collisions T^(0) = 30.4/mb and 
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saturation occurs if the hard cross section is greater than 74 (po/2 GeV) 2 mb. At the LHC, 
gluons alone are sufficient to saturate the transverse area, even with shadowing. For Au 
+ Au collisions at RHIC, the hard cross section must be more than 71 (j9o/2 GeV) 2 mb. 
This condition is not satisfied at RHIC, unless p is lowered to ~ 1 GeV [pTfl . However, 1 
GeV is close to if not within the nonperturbative regime, suggesting that soft physics still 
dominates particle production at RHIC. 

These conclusions depend on the small x behavior of the gluon distribution, the factor 
-fTjet, the cutoff po, and the shadowing parameterization. Transverse saturation does not 
occur at the LHC when the MRST LO set is used if K ]et = 1. An empirical K jet may be 
obtained by comparing model calculations to data, giving some freedom in the value of K- Jet 
for different parton distributions. However, less variation is allowed in the theoretical values 
of -fTjet obtained from the ratio of the NLO and LO cross sections. The theoretical Kj et does, 
however, tends to rise as pr decreases, rendering calculations with p < 2 GeV less reliable. 

Transverse saturation at po = 2 GeV implies that the minijet cross section exceeds the 
inelastic pp cross section, violating unitarity. This is especially a problem for the GRV 94 
LO distributions because of the high gluon density at low x. At very low x then, the proton 
is like a black disk and instead of further splitting to increase the density of partons, the 
partons begin to recombine, acting to lower the density below that without recombination. 
Therefore at very low x, the density of partons should not increase without bound but begin 



to saturate. This recombination corresponds to one picture of shadowing in the proton |16 
A recent HERA measurement of the derivative of the structure function F2 found that at low 
x and Q 2 , dF 2 / 'dlnQ 2 no longer increases, in contrast to the GRV 94 parton densities which 
continue to increase over the range of their validity f2"9 |. The newer MRST distributions 



have been tuned to fit this behavior for Q 2 > 1 GeV 2 . This data implies that the unitarity 
violation in pp interactions is likely an artifact of the free proton parton densities. 

The magnitude of the problem can be gauged by calculating the number of collisions 
suffered by incoming partons. If, on average, a parton collides more than once while crossing 
the nucleus, unitarity violation is a serious problem. The higher the incoming parton x±, 
the more low 2 2 target partons are available for it to interact with, the larger the interaction 
cross section, and the subsequent number of collisions. The minimum x<i depends on p and 
^snn- Since the gluon interaction cross sections are larger than those of quarks, we focus 
on incoming gluons with x = 0.1. The average number of collisions experienced by such 
an initial gluon at the LHC is shown in Fig. ||(a) and (b) for GRV 94 LO and MRST LO 
distributions respectively. The scattering cross section has been multiplied by the nuclear 
profile function T^(6) to give the number of collisions. A gluon can suffer up to an average of 
5 hard scatterings in central collisions with GRV 94 LO and 5=1. It experiences less than 
one collision in the target when b > 5 — 6 fm. Shadowing reduces the severity of the problem 
by decreasing the number of scatterings by ~ 30%. On the other hand, u and u quarks 
with x = 0.1 typically scatter once or less in the target, even without shadowing. With the 
MRST LO distributions, the unitarity violation is less severe, with 1.4 — 2 scatterings per 
central collision for gluons and 0.5 u or u collisions per central event. 

Therefore we might expect that to satisfy unitarity, transverse saturation cannot be used 
as a criteria for determining p and early equilibration by minijet production is unlikely in 
reality. At the lower RHIC energy, unitarity is always satisfied with incoming partons 
experiencing an average of less than one collision. Figure |(c) and (d) shows this for the 
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gluon. The q and q results are considerably smaller. 

The quark rapidity distribution, da q / dy is indicative of baryon stopping due to hard pro- 
cesses. As Tables |I]-|V| and Figs. |§-(7| show, at LHC energies, the GRV 94 LO parton distribu- 
tions predict considerably larger stopping than MRST LO. These homogeneous shadowing 
cross sections can be converted to dN/dy at any impact parameter by multiplying by Tab (6). 
Although both parton distributions predict similar baryon densities at mid-rapidity, GRV 
94 LO predicts about twice as many baryons at large rapidity than MRST LO. Because of 
the unitarity problems and the high gluon density at low x, at LHC the final state baryon 
number, J dy[dqj dy — dqj dy)/3, exceeds the baryon number of the two incoming nuclei. This 
is a clear result of the unitarity violation. Previous works |23| noted this but suggested that 
the problem is reduced if only central rapities are considered, typically \y\ < 0.5. A better 
solution would include a more complete treatment of multiple scattering. However, such a 
calculation involves even more uncertainties. At RHIC energies, the cross sections are lower, 
and baryon conservation is not an issue. The two sets of parton distributions make similar 
predictions, with MRST LO finding a somewhat higher baryon density at mid-rapidity. 

We present calculations covering the entire range of rapidities, even though at the LHC, 
at large rapidity, \y\ > 5, either x\ or x 2 is outside the stated validity range of the parton 
distributions. This range problem could affect calculations at all y since a parton density 
that satisfies the unitarity bound at po will be different at all rapidities since more of the 
low x rise will be subsumed into higher x values to maintain momentum conservation. 

We would like to determine the effects of shadowing on the quantities such as the energy 
density which are important for our understanding of the initial conditions. The initial 
energy density is directly related to the cross section times first Et moment of each flavor, 
(p Q ) (E T ) , which is calculated within a specific acceptance. A crude approximation of the 
acceptance is 

e M ~ { otherwise 

where y max is the highest measurable rapidity. At leading order, the parton pairs are pro- 
duced back-to-back. The Et distribution of each flavor as a function of impact parameter 
is 



2J| 



= I dPT d V2dydzdz' J2 x i F t A ( x uPT,r, z)x 2 F B (x 2 ,p T} b - f,z') 
dE T d 2 bd 2 r 2* J ij= 

1 (da ij ^ kl . 



13 = 
(hi) 



-x (t,u)5(E T -[5 fk e(y) + 8 fl e(y 2 )]p T ) 
1 + Ou { dt 

+^ (u,i)S(E T -l5 fl e(y) + 6 fk e(y 2 )]p T ) ] j . (15) 
Equation (|I~5"D is valid for Et > -Ermin where the ^Tmin required in pp collisions^ is such 



1 A comparison of the LO and NLO jet Et distributions with UA2 data [ 30 1 suggests that below 
Et = 55 GeV, the discrepancy between the calculations and data can be attributed to further 
higher-order corrections or higher- twist effects such as initial and final-state radiation 1 31 ] . 



S 



that a H < a ine i as tic for E Tmin « 2 GeV at 5.5 TeV and 1 GeV at 200 GeV @. Therefore, 
integration over c? 2 r and > E Tmin reduces Eq. flTBl) to the total hard cross section as a 
function of impact parameter 



E 



d<J f (po) 
E T min d 2 bd 2 rdE< 



f dE T 



da(p ) _ H 
2 ^ = a (6 ' Po) • 



(16) 



The last definition in Eq. flIS|) holds for E Tmin — Po, as in Eq. (Q). 

The first moment is obtained by weighting Eq. (|15|) with E T and integrating over 
E T ; we neglect particle masses so that E T = p T , 

^ d^hfr^ = K]et I d PT d V2dydzdz' ]T xiFffap^, f, z)x 2 Ff \x 2 ,p 2 T ,b - r, z') 



(ki) 



1 + S, 



kl 



da ij -+ kl 
di 



(i,u) +6fi 



da ij -* kl 
di 



(u,t) 



(17) 



The first moment is given as a function of rapidity in Figs. |5|-|T2"| both with and without 
impact parameter averaged shadowing for the GRV 94 LO and MRST LO parton densities 
at LHC and RHIC. The average transverse energy given to a particular parton species in a 
central AB collision is then 



E f T (b,p ) 



dgf(p )(El) 
d 2 b 



(18) 



where da^(p ) (E T ) / d 2 b is the integral of Eq. (|T7D over d 2 r. If the nuclear structure functions 



are homogeneous, then the spatial effects factorize and E T (b,p ) is proportional to Tab(P). 
The first Et moment is proportional to the energy density, as we discuss shortly. 
The second moment of each flavor is calculated similarly, 

d0 d^bdh- T ~ = K]ct I d Prdy2dydzdz' ]T xiFf(xi,p T , f, z)x 2 Ff \x 2l p 2 T ,b - f, z) 

(ki) 



x 



p% 



1 + 5, 



ki 



fk 



da ij ^ kl 
di 



(t,u) +S f i 



da ij ^ kl 



+ 



da i3 ^ ff r . da ij ^ ff ' A 
-j {t,u) + -, (u,t) 
dt dt 



dt 

e(y)e(j/2) 



(u,t) 



(19) 



The terms proportional to e(y)e(y 2 ) in Eq. ( |T9|) correspond to only those processes that 
contain identical particles in the final state: qq — > qq, qq — > qq, qq — > gg and gg — > gg. 
These terms are negligible for / = q and q but large for / = g. Indeed, // = gg in Eq. (|19"D 
contributes ~ 30% of the total second Et moment of the gluon. The second moment is 



E 2f (b,p ) 



d 2 b 



(20) 



where da^ \po){E T ) / ' d 2 b is the integral of Eq. fll9|) over d 2 r. For homogeneous structure 
functions, factorization again occurs and E T ^(b,p ) scales with Tab{o)- 
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We now discuss the results characteristic to specific detectors. We will concentrate on 
the coverage around midrapidity, thereby excluding some detector subsystems from our 
consideration. In all cases we assume full azimuthal coverage. At the LHC, there will be 
two detectors taking data with heavy ion beams, CMS [^], optimized for pp studies but with 
a broad rapidity coverage, 7/ max = 2.4, and ALICE |33[, a dedicated heavy ion experiment 
with central rapidity coverage up to y max = 1. We do not include the ALICE forward 
muon spectrometer. The heaviest ions accelerated will be lead. STAR MM and PHENIX 



35] are the two large heavy ion experiments at RHIC, a dedicated heavy ion collider that 



will accelerate ions through gold. STAR has the larger acceptance at central rapidities, 
Umax — 0.9 for the electromagnetic calorimeter, while the central electron arms of PHENIX 
only cover up to y max = 0.350. The PHENIX muon arms will cover more forward rapidities 
but will not increase the coverage at midrapidity except for high mass lepton pairs such 
as those from T decay. The cross sections per nucleon pair and the first and second Et 
moments with and without homogeneous shadowing are given in Tables [I] and [II] in the given 
CMS and ALICE rapidity ranges respectively. At this energy, shadowing can reduce the 
parton yield and the Et moments by up to a factor of two. The corresponding results from 
RHIC are presented in Tables |Hl| and [IV] for STAR and PHENIX. The effect of shadowing is 
much smaller at RHIC than at the LHC. In fact, with S 3 , gluon antishadowing can increase 
the yield relative to S — 1. Recall that the cross sections and moments are all calculated 
with i^ et = 1 and another choice would scale the results correspondingly. 

The effect of the inhomogeneous shadowing is shown for the first Et moment calculated 
with the GRV 94 LO parton densities in Figs. |l| and H for CMS and STAR. The ALICE 
and PHENIX ratios are similar to those shown here. The ratios of the other moments do not 
differ greatly from the first moment. The impact parameter dependence is calculated using 
Eqs. (|3]) and (^). When x lies in the shadowing region, central collisions are more shadowed 
than the average. In the antishadowing region, central collisions are more antishadowed 
than the average. When b ~ Ra, the homogeneous and inhomogeneous shadowing are 
approximately equal, as might be expected from an inspection of Eqs. (|3p and (|j). When 
b ~ 2i?4, the shadowing or antishadowing is significantly reduced. As b further increases, 
the approach to S = 1 is asymptotic. With S p , Eq. ([|), the central shadowing is somewhat 
stronger than with S\vs an d the strength of the shadowing decreases more rapidly when 
b > Ra- At b ~ 2Ra, the ratio with S p is 5% higher than with Sws- A calculation with 
Sr., Eq. (H), would have somewhat stronger shadowing than S p . Since the inhomogeneous 
calculations agree within a few percent, the exact dependence cannot be experimentally 
distinguished; only the presence of inhomogeneity can be detected. This conclusion applies 
to a wide variety of observables M . Because the differences are small, we use only the Sws 
parameterization in the remainder of this work. 

The figures show the ratio of the first Et moment with shadowing included relative to 
S = 1 as a function of impact parameter for q + q, g, and the total, q + q + g. At the 
LHC, quarks and antiquarks are ~ 10% of the total minijet production for the GRV 94 LO 
parton densities and ~ 17% with the MRST LO densities when S = 1. The overall q + q 



2 Since we quote results over full azimuthal coverage, the actual PHENIX cross sections would be 
lower because the central electron arms only cover a fraction of the total azimuth. 
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contribution decreases 1-2% when shadowing is included. At RHIC, the q + q fraction is 
« 19% with the GRV 94 LO densities and « 23% with the MRST LO set. There is again 
only an « 1% variation in the fraction with shadowing included. The ratios in the given 



rapidity intervals with homogeneous shadowing are given by the horizontal lines on Figs. [L3 



and [TJ] and correspond to the ratio of the moments in Figs. ^ and [11] integrated over the 



same rapidity intervals. At the LHC, the ratios are nearly the same for the gluon fraction 
and the total because gluons dominate minijet production. At RHIC energies, since the 
q + q contribution is a larger fraction of the total, the difference between the g ratio and the 
total, shown in Figs. [TBI and [HL is visible, particularly for 5*3 which is shadowed for q + q 



and antishadowed for the gluons. The total remains antishadowed, but less than for gluons 
alone. 

The homogeneous shadowing ratios can also be determined for the zeroth moment, parti- 
cle number, and the second moment of of the Et distribution, from Tables 0-[lV[ The second 
Et moment has a slightly smaller q + q fraction due to the term in Eq. fll9"|) proportional to 
e(y)e(y2) which arises from identical particles in the final state. The dominant contribution 
to this term is gg — > gg, enhancing the overall gluon contribution. 

We now show how the initial conditions for further evolution of the system are impacted 
by shadowing. The initial energy density of each parton species is the ratio of the first Et 
moment to the initial volume 

e{(b,Po) = • (21) 

The total initial energy density is the sum over all species, e, = D/ e f- The initial number 
densities are likewise 

n{ {b , Po) = (22) 

and rii = J2f n {- The energy and number densities are given in Tables |V|- |V11J| , both for 
gluons only and the total minijet yield. Results are shown for both homogeneous and inho- 
mogeneous shadowing at b = where the volume is most clearly defined. Since shadowing 
is stronger in central collisions, the energy density and multiplicity are reduced with S^ s 
relative to the homogeneous case. (See Figs. [13] and [L4] for the impact parameter dependence 
of E T .) At the LHC, inhomogeneous shadowing reduces the energy density by 3-8% at 
6 = 0. At RHIC, the difference is smaller and the energy density may even rise marginally 
at b = with the S3 parameterization. The average energy per particle for a given species is 
e{ /n{ ~ 3 GeV, somewhat larger than p , as can be expected since E T reflects the average 
Pt within the rapidity range. 

These densities can be compared to those obtained for an ideal gas in thermal equilibrium. 
An ideal gas has energy density e t h = 3aT t ^ and entropy density s t h — 4aT t ^ = 3.6n th where 
a = 167r 2 /90 for a gluon gas and a = 47.57r 2 /90 for a three flavor quark-gluon plasmaQ. 



3 N. Hammon et al. also calculated N and Et using spatially homogeneous nuclear structure 
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The initial equilibrium temperatures of such gases are then T th = (e th /3a) 1//4 and the ideal 
energy per particle is 

— = 2.7 T th • (23) 
nth 



We use the results of Tables |Vl- |VIII| with the assumption that = e t h- When 5 = 1, 



T t h ~ 1-07 GeV for gluons only and ~ 840 MeV for a quark-gluon plasma at the LHC with 
the GRV 94 LO distributions and p = 2 GeV. Using the MRST LO results with p = 2 GeV 
yields T th « 860 MeV and 680 MeV respectively. The calculated initial quark-gluon plasma 
temperature is lower than that for gluons because, even though e« is larger for the sum 
of all species, the larger number of available degrees of freedom reduces the temperature. 
Shadowing reduces Th by 10-17% for the gluons and 10-15% for the total with the largest 
effect due to Si and the smallest from S3 with its antishadowing. At RHIC, the equivalent 
temperatures extracted with po = 2 GeV are smaller, 410 MeV for gluons and 330 MeV for 
a quark-gluon plasma with S — 1. The reduction due to shadowing is 5% or less — in fact 
a slight enhancement is possible because of the antishadowing in S3. The temperatures are 
virtually independent of the parton distributions at this energy since the two sets are very 
similar in the x range of RHIC. The temperatures estimated for RHIC are lower than those 
obtained elsewhere. This difference will be discussed in the next section. 

These equivalent equilibrium temperatures are only approximate because they depend 
on the rapidity range over which Th is calculated. The extracted temperature rises as the 
rapidity range decreases because the antiquark and gluon distributions are maximal at y = 0. 
The fact that the width of the slice affects Th shows that thermalization in the collision 
is incomplete. To study this further, we can compare these results with expectations from 
the ideal gas. The GRV 94 LO gluon temperature, T t h ~ 1-07 GeV, satisfies Eq. (|23|) when 
S = 1, i.e. 



nf(6 = 0,po) 



«-fL = 2.7T th . (24) 
5=1 < h 



This equation suggests that, even if a quark-gluon plasma is far from equilibrium, the gluons 
might equilibrate quickly, around n ~ 0.1 fm, even without the secondary collisions required 
for isotropization. However, even this suggestion only holds at LHC energies without shad- 
owing. Shadowing drives the result away from equilibrium so that 



4( b ,Po) _ 4h 



nf{b,p ) 



>^L = 2 .7Tt h . (25) 



Note however that taking K- }et = 1.5 increases all the extracted temperatures by « 10%, 
bringing the shadowed results closer to the ideal in Eq. fl23"|) . Reducing p for S 7^ 1 at the 
LHC or for any scenario at RHIC would also increase the extracted temperatures so that the 



functions at RHIC and LHC |36|. Since they take i£j e t > 1, they find larger energy densities and 
effective temperatures than we do here. They also neglected the unitarity problem in their LHC 
estimates. 
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gluon result would again appear to equilibrate. The quarks alone or the q + q + g total will 
not come to equilibrium, even when p Q is reduced, due to the lower equivalent temperature. 
In any case, po cannot be set arbitrarily low for perturbative QCD to be valid. In addition, 
p should not be a strong function of energy and should be independent of the shadowing 
parameterization. Shadowing thus reduces the likelihood of fast thermalization, even at the 
LHC where the conditions are most favorable. 

Given these uncertainties, one can nevertheless obtain an approximate lower bound on 
the produced particle multiplicity. In an ideal longitudinally expanding plasma, the energy 
density evolves following |J7] 

£ + l±^ = (26) 

GST T 

where P is the pressure and r is the proper time. There are two extreme solutions: free 
streaming, with P = 0, leading to e ~ r _1 , and ideal hydrodynamics, P = e/3, where e ~ 
r -4/3 rj^ e j ower jj m jt Q f multiplicity is obtained from ideal hydrodynamics where the system 
is treated as it it were in thermal equilibrium at 73 = l/po = 0.1 fm and expands adiabatically 
with t. Then the initial entropy determines the final-state multiplicity, neglecting final-state 
interactions, fragmentation, and hadronization. If only minijet production contributes to 
the final-state multiplicity, the total number of particles in a specific detector's central 
acceptance is then |38] 



dN 4 



dy 3.6 



Ti-KR\a I E T (\y\ < 

1/max ) 



3, 



(27) 



27 [ 2|/ max J 
Equation fl2"7|) suffers from some uncertainty due to the t*' 4 ~ p l Q A dependence in the 



volume besides the dependence on p in Et- In a complete calculation, the variation with 
Po would be compensated by a corresponding variation in the soft component, as discussed 
later. With the GRV 94 LO distributions at the LHC the total dN/dy at y = from 
minijets is « 4000 - 6000 or about 2700 - 4000 charged particles, ^2/3 the total dN/dy. 
Shadowing reduces the number of charged particles to ~ 1800 — 2600. With the MRST LO 
distributions, the total dN/dy\ y=0 « 2000 — 3500 without shadowing and 1400 — 2600 with 
shadowing. With inhomogeneous shadowing, the LHC multiplicity drops 2-5% for collisions 
at b = 0. The gluon Et moment dominates the total and drives the rapidity distribution, 



as can be inferred from Figs. |9] and [Ty. We find total minijet multiplicities of 220 — 350 
without shadowing and 200 — 360 with shadowing. The larger dN/ dy with shadowing is a 
result of the antishadowing in S3. Since soft production is large at RHIC, the total dN/dy 



found here is considerably lower than predicted by some event generators 39 



IV. CORRELATION BETWEEN E T AND IMPACT PARAMETER 

Thus far, we have discussed the dependence of shadowing on the impact parameter, a 
quantity which cannot be directly determined in a heavy ion collision. However, although 
the impact parameter is not measurable it can be related to direct observables such as the 
transverse energy, Et 0,^4| ■ The transverse energy is summed over all detected particles in 



the event with masses mk and transverse momenta pTk so that Et = 'EkJml + p Tk . Besides 
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an Et measurement, it is also possible to infer the impact parameter by a measurement 
of the nuclear breakup since the beam remnants deposited in a zero degree calorimeter are 
also correlated with the impact parameter [TO. A measure of the total charged particle 



multiplicity, proportional to Et [0 , could also refine the impact parameter determination. 

The transverse energy contains "soft" and "hard" components. The "hard" components, 
calculated in the previous section, arise from quark and gluon interactions above the cutoff 
Po = 2 GeV. "Soft" processes with pt < Po are not perturbatively calculable yet they can 
contribute a substantial fraction of the measured Ej- at high energies (and essentially the 
entire Et at CERN SPS energies). These processes must be modeled phenomenologically. 
Our calculation of the total Et distribution follows Ref . |2|J . We assume that the soft cross 
section, a^ ' , is equal to a™ elastlc , the inelastic pp scattering cross section. The hard part of 



the Et distribution can be expressed as 



da 



H 



dEn 



N\ 



N 



exp 



N 

Y[dE T i 



1 da, 



H 

pp 



, a^ p dE Ti 



5(E T — S i=1 .ETi 



(28) 



The average number of hard parton-parton collisions is defined in Eq. 



2Ra, N is large and da H jdET can be approximated by the Gaussian p 



da 



H 



dEn 



d 2 b 



27ia 2 E H 



(b) 



exp 



T 



E H T {b)f 



2o%" 



For most b < 



(29) 



where the mean Et, E t (6), is proportional to the first moment of the hard cross section, 

E"(b) = Z f E f T (b, Po ). (30) 
The standard deviation, a^(b) is computed from the first and second moments, 



4>) = E4 7 (Mo)-4tA 

T a (b,p ) 



(31) 



see Eqs. ( PD and (|2"0"|). The impact parameter averaged values of the hard cross section 
and its first and second Et moments correspond to the "total" values in Tables ||-[IV] for the 
specified rapidity coverages of the four detectors. Note that these moments are a lower bound 
on particle production from hard processes because hadronization has not been included. 

The soft component is usually taken to be proportional to the number of nucleon-nucleon 
collisions, 



N\b)=T AB {b)a 



s 

pp ' 



(32) 



where a^ p 
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Since a pp 



40 mb at RHIC and may increase to 60 mb at the LHC 
depends only weakly on the collision energy, the hard and soft components are assumed to 
be separable on the pp level and thus independent of each other at fixed b [24]. The soft 



component may be computed using the first moment, E T , and second moment, E\ , of the 
soft E T distributions, obtainable from lower energy data [Q. At the SPS, ^/snn = 19.4 
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GeV, o. 



vv 



32mb, En 



15 mb GeV and 



50 mb GeV 2 ]2J for \y\ < 0.5. We 

are independent of impact parameter and scale with energy as 
<jp V and linearly with rapidity acceptance. The resulting first and second moments for the 

could be scaled by the 
21. However, at these 



g c 

assume that E T and E T 



four detectors are given in Table |X|. Alternatively, E T and E\ 
charged particle production rate in the selected rapidity interval 
higher energies, the charged particle distributions will have a strong contribution from hard 
production which could lead to double counting of the total rate. If the SPS multiplicity 
distribution is used, then the effect of the rising cross section will dominate. 

The total Et distribution is a convolution of the hard and soft components with mean 
and standard deviation 



E T (b) 



E*(b)+T A 



- / - J-AB(b)E T 
a E 



\b)+T AB (b)a 2 E 



(33) 
(34) 



The standard deviation for the soft component, a E , is 



^-2 s 



■pS 2 



a 



pp 



(35) 



We do not assume that the second moment E T 



is equivalent to the standard deviation as 



in some previous calculations 



Some caveats related to the soft Et contribution should be mentioned. The second mo- 
ment is system dependent^, perhaps because the fluctuations are concentrated in the central 



region, making E T sensitive to the acceptance |[44|| . There also may be some contamination 



from hard processes. Additionally, soft processes may also be subject to a form of shadowing 
due to large mass diffraction |45| , analogous to the multiple-scattering picture of shadowing 
except that it affects soft interactions. If correct, the soft component would also be reduced 
and the soft and hard interactions would have a similar impact parameter dependence. Thus 
the soft component is only accurate to the 20% level at best. 

At the LHC, the hard component is an order of magnitude larger than the soft part. 
This can be seen from a comparison of the homogeneous shadowing first and second Et 
moments in Tables |-[IV| with E T and E\ in Table [IX]. The results are directly comparable 
because the first and second Et moments in all the tables are given per nucleon pair. At the 

5 

LHC, the moments from minijet production are 6 — 10 times larger than E T with the GRV 

__ g 

94 LO parton densities and 3 — 7 times larger than E T when calculated with the MRST LO 
parton densities. Total particle production is then dominated by minijet production. With 
soft production included, the estimated dN/dy in Sec. Ill would be increased by 5 — 22%, 
less than the change due to shadowing. 

At RHIC however, E T is 1.3 — 2.2 times larger than the first Et moment, depending 
on the parton densities and shadowing parameterization. Thus the soft contribution to 
the total Et is still somewhat larger than the hard contribution. When soft production 



4 See oj in Table 8 of Ref. [43|. In lighter targets to is significantly different than in heavy targets. 
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Q £J _____ 

is included in the estimated dN/dy by adding TAB(b)E T to E T (b,p ) in Eq. (j27|) , dN/dy 
at b = could increase by a factor of 1.9 — 2.4, up to 680 — 750 particles. Likewise, the 
extracted initial temperature assuming thermal equilibrium would be ~ 20% higher when 
the soft contribution is included and could reach « 500 MeV with £7=1, consistent with 
previous predictions f46|. If soft production is also affected by shadowing f45|, then the soft 



contribution to RHIC central collisions would be reduced and the hard and soft components 
would be more in balance. 

The Ex distributions, with homogeneous and inhomogeneous shadowing are shown for 
each detector in Figs. p~5|4T8| - The hard component is calculated with the MRST LO distribu- 
tions. In each case, we show the change in the Et distribution due to shadowing in the most 
central collisions, b < 0.2Ra, semi-central collisions, 0.9Ra < b < I.IRa, and the entire b 
range. The maximum Et is reduced 30 — 40% at the LHC because the hard component, 
Eq. fl33|) , dominates the average Et- At intermediate impact parameters, the Gaussian, 
Eq. (|29"D , is narrowed by shadowing. At RHIC, since the hard and soft components are 
comparable, the maximum Et is shifted by only ~ 7% when shadowing is included. Indeed, 
for £7 3 , since shadowing enhances the Et moments of the hard component, the maximum 
Et is slightly increased. If the GRV 94 LO distributions are used in the calculation of the 
hard part, the total Et at the LHC is nearly twice as large and the shadowing effects are 
stronger. The RHIC results are essentially unaffected by the choice of parton distribution 
since the Et moments do not depend strongly on the parton distribution, see Tables |I| and 



IV. 

These results depend on Kj et since the hard E T is proportional to Kj et . At the LHC, 
Et scales nearly linearly with K )et since hard interactions dominate there. At RHIC, the 
increase would be smaller, since only 30%-50% of the Et comes from hard processes; if 
Kj et = 1.5, then the maximum Et rises by 20%. Similar results were found in Ref. [pfifl . 

The change in the Et distribution due to shadowing is not equivalent to scaling Et by 
a constant. The shape of the distribution is also modified because central and peripheral 
collisions are affected differently. The shape change is small at RHIC, but clearly visible for 
the LHC. Figures [15] and [16] show that the shadowed distributions are enhanced over £7=1 
for Et ~ 6 TeV and 3 TeV for CMS and ALICE respectively. If soft production is also 
affected by shadowing |[45|| , the shape change may be larger for RHIC. 

For semi-central through central collisions, the transverse energy-impact parameter cor- 
relation is relatively easy to determine but in very peripheral collisions, the entire transverse 
energy could arise from a single hard collision which produces e.g. a J/ip or a Drell-Yan 
pair. Then, the simple Gaussian approximation to Eq. would break down. 



V. DRELL-YAN, J/ip AND T PRODUCTION 



We now study the effect of inhomogeneous shadowing on the production of hard probes. 
As examples, we consider Drell-Yan and quarkonium production. We have previously studied 
the production of charm and bottom quarks at these energies ||. We have also considered 
shadowing effects on J/ip and Drell-Yan production at the SPS, as well as their ratio as a 
function of Et [57]. However, at the SPS, Et is dominated by the soft component and is 
proportional to the number of participants |4T |. We do not include final-state absorption 
effects on quarkonium production. 
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These calculations are done at leading order to be consistent with our calculations of 
minijet production. The LO cross section for nuclei A and B colliding at impact parameter 
b and producing a vector particle V (quarkonium or 7*) with mass m at scale Q is 

da v ^ r , . ,„„, „ ^ , „ D , „o t* ^ ,, do. 



r -* da v - 

^ J dzdz>F A (x 1 ,Q^r,z)Ff(x 2 ,Q\b-r,z')^-^ , (36) 



dydm 2 d 2 bd 2 r ~f J ' J dydm 

where <r^ is the partonic — > V cross section and the parton distributions are defined in 
Eq. ©. 

The LO Drell-Yan cross section per nucleon must include the nuclear isospin since, in 
general, + a° Y ^ a° Y + a° Y , 

fc«^«£ = «,^ (37) 

q=u,d,s 

where Z4 and Na are the number of protons and neutrons in the nucleus. We assume 
charge symmetry, = /| = /" etc., in the nuclear environment. In Eq. ([37]), Xi j2 = 
Qe ±y / ^snn and Q = m. The factor -Kexp, typically 1.7 — 2 for fixed-target Drell-Yan 
production, accounts for the difference in magnitude between the calculations and the data. 



Figures [19| and |2(] show the influence of shadowing on the Drell-Yan mass distribution, 
calculated with the MRST LO parton distributions. The ratio of the inhomogeneously shad- 
owed mass distribution to that for 5* = 1 are shown in several impact parameter bins, along 
with the homogeneous shadowing ratios, in the rapidity coverage given for the ALICE and 
PHENIX central detectors. The corresponding ratios for CMS and STAR are quite simi- 
lar. Ratios are presented for the most central collisions, b < 0.2Ra, semi-central collisions, 
0.9Ra < b < I.IRa, and peripheral collisions, 1.9Ra < b < 2.1Ra- In the most central 
collisions, the inhomogeneous shadowing, with Eq. ([3]), is somewhat stronger while in the 
most peripheral collisions, it is much weaker. In each case, the S3 parameterization gives 
the smallest effect. At the LHC, evolution is also most apparent with this parameteriza- 
tion. A shortcoming of the limited Q 2 evolution of the 5*2 parameterization is obvious in 
Fig. [TD| — the evolution is evident up to m = 10 GeV after which the 10 GeV values of the 
valence quark, sea quark, and gluon shadowing ratios are used at all higher masses. Above 
10 GeV, the ratios with the Si and 62 parameterizations are then similar. The Si results 
change very slowly with mass because they lack Q 2 evolution. At the lower RHIC energy, 
the 10 GeV Q 2 cutoff in S2 is less obvious because the x values are larger, in a region where 
shadowing is small. At RHIC, shadowing of the most peripheral collisions predominantly 
occurs for masses below 8 GeV. At this energy the largest mass pairs are antishadowed. The 



antishadowing is weakened in peripheral collisions, see Fig. E0 



Since the next-to-leading order, NLO, Drell-Yan cross section includes Compton scat- 



tering with an initial gluon P5| , it is possible that shadowing could change significantly at 
NLO, especially with the S% and S3 parameterizations. We have therefore also calculated 
the Drell-Yan cross sections at NLO with all the homogeneous shadowing parameterizations 
and found that the ratios do not change significantly when the NLO terms are added. There 
is a 3-4% difference in the ratios with shadowing at LO and NLO in Pb+Pb collisions at 5.5 
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TeV and 0.5-1% in Au+Au collisions at 200 GeV. This should not be too surprising since 
the theoretical K factor is small, K t ^ = ct^lo/^lo ~ 1-2 a t RHIC and 1.1 at the LHC. 
The effect of shadowing on the higher order contributions must then be less than K t ^, small 



compared to the uncertainties in the shadowing model, as can be seen from Fig. £L 



Figures |22] and ^3] show the rapidity dependence of the shadowing for Drell-Yan pro- 
duction when 4 < m < 9 GeV. The homogeneous and inhomogeneous results are again 
compared in central, semi-central, and peripheral collisions. The S 2 parameterization pro- 
duces the strongest shadowing because the sea quark ratio is lower at small x than Si and 
5*3, see Fig. Q All the LHC ratios increase with rapidity because x 2 remains small while X\ 
increases to ~ 0.1 at y ~ 5. Recall that around x\ ~ 0.1 Si shows antishadowing, S 2 ~ 1 
for sea quarks, and the sea quark distributions are shadowed with the S3 parameterization. 
Thus the change in the shadowing ratios as a function of y is smallest with S3. As y and 
X\ increase, the shadowing, antishadowing, and EMC regions are traced out. However, at 
forward rapidities, x 2 < 10~ 4 so that the cross section ratios are always significantly less 
than unity. 

At RHIC, the ratios decrease with rapidity. Both x\ and x 2 are in a region where all the 
parton densities are shadowed at y — but, as the rapidity increases, x 2 decreases to the low 
x saturation region while the X\ values enter the EMC region. The resulting convolution is 
then lower at large y than at central rapidities. Since the Drell-Yan cross section is calculated 
in the interval 4 < m < 9 GeV, some influence of Fermi motion is apparent at the largest 
rapidities because x\ ~ 0.9 when m = 9 GeV and y = 3. 



The effect of the inhomogeneity is shown more fully in Figs. [24| and We have chosen 
two different mass ranges, 4 < m < 9 GeV and 11 < m < 20 GeV, between the J/ip 
and T resonances and above the T family respectively. The similarities between the CMS 
and ALICE predictions at the LHC and the STAR and PHENIX expectations at RHIC 
are obvious in these figures. In the range 4 < in < 9 GeV, shadowing is expected at all 
masses. In the larger mass region, the similarity between the Si and S2 parameterizations 
above 10 GeV are visible in the CMS and ALICE plots. For completeness, the LO Drell-Yan 
production cross sections per nucleon pair for both mass ranges are shown in Table [X] with 
and without homogeneous shadowing. Recall that the theoretical K factor between the LO 
and NLO cross sections is ~ 1.2. 

We now consider shadowing in J/ip and T production using two models that have been 
successfully employed to describe quarkonium hadroproduction. The first, the color evap- 
oration model, treats all quarkonium production identically to QQ production below the 
MM threshold, where M represents the lightest meson containing a single heavy quark Q, 
neglecting the color and spin of the produced QQ pair. The non-relativistic QCD approach 
expands quarkonium production in powers of v , the relative Q-Q velocity within the bound 
state. In this model, the produced QQ pair retains the information on its color, spin and 
total angular momentum, requiring more parameters than the color evaporation model. 



In the color evaporation model [49 



/,"(*!, Q 2 )ff K Q 2 )^H = Fc,sK th (/"(a*, Q 2 )L N (x 2 , Q 2 ) 



s 

2\ " 



E [f^^Q^f^^Q^ + f-^x^Q^f^x^Q 2 )]^^} , (38) 



q=u,d,s 
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where C and B represent the produced charmonium and bottomonium states. The LO 
partonic QQ cross sections are defined in [[[(J and s = X1X2S. The fraction of QQ pairs 
below the MM threshold that become the final quarkonium state, Fqb is fixed at NLO |19| . 
The factor matches the LO cross section to the NLO result. Together, the multiplicative 
factors Fc,b and K t ^ reproduce the pp data in magnitude and shape. For J /if) production, 
we use m c = 1.3 GeV and Q = m c with the GRV 94 LO distributions and m c = 1.2 GeV 
and Q = 2m c with the MRST LO densities [J49|. For T production, we take m& = Q = 4.75 



GeV with both sets of parton distributions. 

The J /if) cross section ratios in the color evaporation model are given as a function of 
rapidity at LHC and RHIC in Figs. [26] and |27| respectively. At both energies, the S\ and 
S 2 results are very similar because the product of the S\ shadowing ratios and the S 2 gluon 
shadowing ratios at Q = 2m c = 2.4 GeV differ by only 1 — 2% over a wide range, 5 units 
of rapidity at the LHC and 2.5 units at RHIC. The ratios with the S3 parameterization 
are larger than with the Si and S 2 parameterizations. This is due to the nature of the S3 
parameterization: at low x and \y\ there is less gluon shadowing and at large x and \y\ the 
gluon antishadowing is stronger than in Si and S 2 . These effects are also obvious in the 
rapidity- integrated impact parameter dependence shown in Fig. ^S. 

The J /if) results in the color evaporation model are rather sensitive to the choice of parton 
distributions. This sensitivity arises from the rather low m c compared to the initial scale of 
many parton distributions. The initial scale of the MRST LO densities is Qq ~ m c w 1.1 GeV 
suggesting Q = 2m c is an appropriate choice. Because the initial scale in the GRV 94 LO 
densities is Qq ~ m c /2 ps 0.63 GeV, we use Q = m c . Choosing the scale proportional to m c 
is somewhat more consistent with the calculations of the Drell-Yan and minijet production 
cross sections. However, the light charm quark mass precludes this choice for the MRST LO 
densities. We have displayed the results with the MRST LO densities. If the GRV 94 LO 
densities are used, the shadowing is somewhat stronger at both energies and the Si and S2 
results are different. 

Figures. [29| and [30] show the shadowed T cross sections, relative to S = 1, as a function 
of rapidity in several impact parameter regions. The Si and S 2 parameterizations now 
differ due to the evolution of the S2 parameterization. The Si parameterization, without 
evolution, gives an T ratio only slightly different from that of the J /if) at y — for the LHC 
energy because as X2 changes from 4.4 x 10~ 4 for the J /if) to 1.7 x 10~ 3 for the T at y = 0, 
Si is nearly constant, see Fig. ||. The peak at y ~ 4.3 with Si appears as xi goes through 
the antishadowing region to the EMC region. While the maximum in the shadowing ratios 
occurs at similar rapidities in J /if) production, y ~ 5 for S 3 and y ~ 5.5 for Si and S 2 , 
the T ratios peak at y ~ 3.5 for S 3 , 4.5 for Si and y ~ 5 for S 2 . In fact, now the S 2 
and S3 ratios are similar at the LHC. The larger gluon antishadowing associated with J /if) 
production is reduced at the larger bottom mass. At RHIC T shadowing is further reduced 
relative to the J /if} than at LHC. In contrast to Fig. |27], the ratio decreases with increasing 
y over all rapidity. Note also that T production is restricted to a narrower range than the 
J/if> because the T is heavier. Little T shadowing is observed with S 2 while S3 exhibits 
strong antishadowing at y = since xi = x 2 = 0.048. The T results are less dependent on 
the choice of parton distributions than the J /if). This set of parton distributions is weaker 
than that of the J /if). This is because m& > Qq in both sets so that we choose m& = Q, 
eliminating the ambiguity in scale due to the small charm quark mass in J /if) production. 
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The impact parameter dependence of T production is shown for the central rapidity 
coverages of the LHC and RHIC detectors in Fig. |3T]. These T ratios are much more 
dependent on rapidity than the corresponding J/ip ratios. Since the largest shadowing or 
antishadowing occurs in the central region, a stronger relative ^-integrated effect is observed 
in the detectors with the narrowest rapidity acceptances. This is particularly obvious for 
the S3 parameterization in PHENIX with respect to STAR. 

The effects of shadowing on quarkonium production in the color evaporation model are 
unchanged between LO and NLO ||51|| . Even though at NLO quark-gluon scattering also 
contributes to quarkonium production, the fraction of the total production cross section due 
to this new channel is not large enough at these energies to change the shadowing effects. 

The non-relativistic QCD, NRQCD, approach is an extension of the color singlet model 
| 52|| which requires J/^'s to be produced with the correct color and total angular momentum. 
The color singlet model predicts that high J/ip production occurs dominantly through 
Xcj decays because direct J / ip production required a hard gluon emission on a perturbative 
timescale. The NRQCD model does not restrict the angular momentum or color of the 
quarkonium state to the lowest allowed color singlet state. Then, e.g. a J/ip may produced 
as a 3 Po color octet which hadronizes through the emission of nonperturbative soft gluons. 

The rapidity distribution of the final-state C or B is 

f?{^Q 2 )f?{^Q 2 ) d ^ B 



dy 

= £E/^2<% - \\n (^))// v (xi,g 2 )/i v (x 2 ,g 2 )c^ [n] (^) . (39) 

The sum over i and j includes up, down, and strange quarks and antiquarks as well as gluons 
since in NRQCD e.g. the process (q + q)g — > Xc\X also contributes to J/ip production. The 
expansion coefficients C % ^,. are calculated perturbatively in powers of a s (Q 2 ) up to a 3 and 

the nonperturbative parameters (0^ ,B ) describe the hadronization of the quarkonium state. 
The expressions for the cross sections and the values of the nonperturbative parameters can 
be found in Ref. fl55| . Since {0%' B ) were fixed using the CTEQ 3L parton densities |54| with 
m c = 1.5 GeV, rrib = 4.9 GeV, and Q = 2mQ, we use this set with the same uiq and Q 
values to be consistent with fixed target cross sections 



The total J/ip cross section includes radiative decays of the Xcj states and hadronic 
decays of the tp', 

^ = ^ + £b^ mX)^ + B& W | . (40, 

Likewise, the total T cross section includes radiative decays from x fe j(lP) and x 6 j(2P) 
states and hadronic decays from the T(2S") and T(3S) states. We have not included ra- 
diative decays from the proposed xw(3-P) states since their branching ratios to the lower 
bottomonium states are unknown. Then 

dar - daf + E B( XbJ (iP) - rx)%^ + b cS (?(2S) - Tx) dari2S) 



dy d y j= dy dy 

f E B eS ( X u(2P) - TX)^f^ + B cS (T(3S) - TX)^£) . 
J^o dy dy 
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Note that as well as the direct decays of the higher bottomonium states to the T, a final- 
state T can be produced by a chain of hadronic and radiative decays. In the case of e.g. 
the T(3S f ), decays to T(2S) and T are of the same order as decays to the Xbj(2P) states. 
The branching ratios above the Xbji^P) states are labeled as B e Q to indicate that direct as 
well as chain decays are included in the total branching ratio. The perturbative part of the 
production, Cggr^ is the same for the T, T(2S), and T(35) states and for the Xw(l-P) and 
Xbj{2P) states. Only the parameters (O^) change. The complex feeddown of the higher 
bottomonium states to the T requires more parameters than J/i/j production. 

In contrast, in the color evaporation model, the rapidity distributions of all states are 
assumed to be the same, thus e.g. Fj/^ in Eq. includes the Xcj and ip' decay contributions 
given explicitly in Eq. (fX0|). 

Two differences between the NRQCD and color evaporation approaches are relevant here. 
The first concerns the x values probed. Since the color evaporation model integrates over 
QQ pair mass up to the MM threshold, it averages over the x range 2mg/ \/snn < x < 
2rriM I y/SNN- The pair mass integration also includes limited Q 2 evolution in the parton 
densities and the shadowing parameterizations. The NRQCD formulation selects specific 
X\ and £2 values for some of the states and only involves a convolution over x for color 
singlet production of e.g. gg — > J/ip,Xci,Xc2 and g(q + q) — > Xa- Additionally, production 
is at fixed Q 2 for all states. The second difference is the g(q + q) contribution to NRQCD 
production, absent in the color evaporation model. 

We show NRQCD results for J/ip production in Figs. and [33|. Since the Si parame- 
terization is flavor and Q 2 independent, these results are least influenced by the production 
model. The differences between the models are most obvious at RHIC where the qq contri- 
bution is ~ 5% of the color evaporation cross section and ~ 1% of the NRQCD cross section. 
The g(q + q) contribution is ~ 3 — 4% of the NRQCD cross section. Since the gluon is anti- 
shadowed at RHIC, significantly less shadowing can be expected in the NRQCD model than 
in the color evaporation model. The relative reduction in shadowing is particularly obvious 
for the S3 parameterization in Fig. |33] where the cross section ratio is ~ 0.95 over 1.5 units 
of rapidity where x is antishadowed. At larger rapidity, x is in the EMC region and the S3 
gluon ratio decreases again, as shown in Fig. |2|. The S2 ratio is generally flatter because the 
gluon ratio is not reduced in the EMC region. The difference between the two approaches is 
significantly smaller at the LHC where the qq contribution is less than 1% for both models 
and therefore plays practically no role. 

The impact parameter dependence of shadowing in the NRQCD approach on J/ip pro- 



duction is shown in Fig. |34|. The difference between shadowing in this model and in the 
color evaporation model seen in the rapidity distributions is obvious here as well. 

The effect of shadowing on T production in the NRQCD approach is shown on the 
rapidity distributions in Figs. ^ and [36| and on the impact parameter dependence in Fig. |37 . 
The same trends seen in the color evaporation model are observed here except that shadowing 
or antishadowing effects are reduced for NRQCD production. Here, the larger h quark mass, 
4.9 GeV, and scale, Q = 2m&, reduce the magnitude of the shadowing. The importance of 
qq annihilation in the color evaporation model relative to the qq and g(q + q) contributions 
in NRQCD affects the shadowing. The g(q + q) component in NRQCD is 1% or less of the T 
cross section at both RHIC and LHC. The higher quark mass probes larger x values where 
the qq contribution is larger. At RHIC, qq contributes 13 — 16% of the total T cross section 
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in the color evaporation model compared to 36 — 56% of the total T cross section in the 
NRQCD approach. The larger fraction of T production by qq annihilation in NRQCD is 
due to the large octet XbJ contribution. 

The integrated J /if) cross sections per nucleon pair for both models are shown in Table [XI] . 
The factor K t ^ is included for the color evaporation model while the NRQCD parameters 
are fit to the measured cross sections at LO. The 5* = 1 cross sections agree within 5 — 
7% at RHIC and within 15% at the LHC. The NRQCD results are lower than the color 
evaporation results at RHIC but the NRQCD cross section grows faster with energy than 
the color evaporation cross section. This behavior can be attributed to the different small x 
behavior of the MRST LO and CTEQ 3L parton densities. With homogeneous shadowing, 
the differences are more striking, as reflected in Figs. |26|-[34, 



Table [Xil| shows the integrated T production cross sections per nucleon pair for both 



models. The theoretical K factor is included for the color evaporation model [fj9| . The 
NRQCD parameters have been fit to fixed target T production data. The two model T 
cross sections do not agree as well as do those of the J/ip. Reasons for this disagreement 
might include the greater number of T parameters needed to fit a more limited set of data 
or the absence of the possible Xftj(3-P) decays in this calculation. 

Finally, we mention one caveat concerning quarkonium production. Since the initial 
quarkonium state is typically a color octet and obtains its final-state identity in a later 
soft interaction, it is conceivable that production and conversion occur far enough apart 
in position space for the strength of the apparent shadowing to be different. However, if 
shadowing is considered to only affect quarkonium at the production point, this separation 
is insignificant. In any case, this separation is a much bigger issue in pp interactions, where 
the two points must be quite close. 



VI. DISCUSSION AND CONCLUSIONS 



We have studied the effect of shadowing and its position dependence on particle pro- 
duction in nucleus-nucleus collisions at RHIC and LHC energies. Shadowing can reduce 
the minijet yields by up to a factor of two at the LHC. Assuming that hard production 
dominates the determination of the initial conditions and that the high minijet yield leads 
to equilibration, the initial energy density and apparent temperature can be significantly re- 
duced. Fast equilibration is unlikely, even for the gluons alone, when shadowing is included. 
The change in the initial conditions due to shadowing is considerably smaller at RHIC, on 
the order of a few percent, less than the change in the initial conditions when soft production 
is included. We have compared the initial conditions in central collisions with homogeneous 
and inhomogeneous shadowing and found the difference to be small. The inhomogeneity of 
the shadowing becomes more important in peripheral collisions. We have also showed the 
shadowing effects on the E? distributions for the central rapidity acceptance of the major 
detectors at the LHC and RHIC. We note that our results at RHIC are more stable with 
respect to changes in the parton densities than at the LHC where the small x behavior 
of the gluons can lead to unitarity violations, the size of which depends strongly on the 
chosen parton densities. Since we assume Kj et = 1, we have been very conservative in our 
estimates of the initial conditions. Nonetheless, once unitarity is satisfied at the LHC, the 
hard component is likely be reduced judging from the difference between the GRV 94 LO 
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and the MRST LO cross sections. Thus lower number and energy densities may be expected 
for LHC collisions. 

Finally we have studied the effects on the J/ip, T, and Drell-Yan yields. A careful 
measurement of the J/tp, T, and Drell-Yan rates as a function of rapidity can help distinguish 
between shadowing models as well as the quarkonium production mechanism since the color 
evaporation and NRQCD approaches lead to quite different shadowing patterns. Because 
there is typically a larger shadowing effect on quarkonium production in the color evaporation 
model than on Drell-Yan production, e.g. the J/ip to Drell-Yan ratio would be smaller than 
that expected for S — 1. On the other hand, the NRQCD approach predicts the reverse- 
the J/ip to Drell-Yan ratio may be larger than expected when S = 1. Since the effect of 
shadowing depends on the Drell-Yan pair mass, if the Drell-Yan yield is to be used as a 
baseline to compare the yield of other hard probes, the rates should be measured directly 
in the mass region of interest rather than relying on calculations to extrapolate into an 
unmeasured region. 

One key test of the impact parameter dependence of shadowing is the slope of the Drell- 
Yan mass distribution; if shadowing varies with position, the slope of the distribution should 
depend on Et- If the slopes are significantly different for central, intermediate and peripheral 
collisions, this would be a clear demonstration that shadowing depends on position. The 
only complication may be due to parton energy loss before the hard interaction. 
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TABLES 
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S3 
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TABLE I. The minijet cross section, Eq. @, first and second moments of the transverse energy 
distribution, Eqs. ( |T7| ) and ( |l9| ) respectively with po = 2 GeV, integrated over b and r and divided 
by AB, within CMS, \y\ < 2.4. Results for both sets of parton distributions used are separated 
into contributions from quarks, antiquaries and gluons as well as the total. The calculations are 
done without shadowing, S = 1, and with shadowing parameterizations Si, S2, and S3. 
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TABLE II. The minijet cross section, Eq. (BI), first and second moments of the transverse 



energy distribution, Eqs. (17) and ( |19|) respectively with po = 2 GeV, integrated over b and r and 
divided by AB, within ALICE, \y\ < 1. Results for both sets of parton distributions are separated 
into contributions from quarks, antiquaries and gluons as well as the total. The calculations are 
done without shadowing, S = 1, and with shadowing parameterizations Si, S2, and S3. 
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TABLE III. The minijet cross section, Eq. (|6|), first and second moments of the transverse 
energy distribution, Eqs. (|l7|) and (|l9|) respectively for p$ = 2 GeV, integrated over b and r and 
divided by AB, within STAR, \y\ < 0.9. Results for both sets of parton distributions are separated 
into contributions from quarks, antiquaries and gluons as well as the total. The calculations are 
done without shadowing, S = 1, and with shadowing parameterizations Si, S2, and S3. 
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GRV 94 LO 






MRST LO 




c 
o 


q 


Q 


9 


total 


Q 


q 


9 


total 










2a f (po) (mb) 










1 


n on 


0.18 


2.28 


2.75 


0.25 


0.16 


1.81 


1. 11 




(J. 24 


0.15 


1.86 


2.25 


0.20 


0.13 


1.48 


1.81 


bo 


U.z4 


0.14 


1.91 


2.29 


0.21 


0.13 


1.53 


l.o ( 


<->3 


9R 


0.17 


2.32 


2.77 


0.24 


0.15 


1.85 


9 9/1 








a 


% )(i4> ( m b GeV) 








1 


n oi 


0.50 


6.01 


7.32 


0.69 


0.44 


4.85 


K HQ 

o.yo 


c 

bl 


U.OO 


0.41 


4.98 


6.07 


0.58 


0.37 


4.04 


A on 


C 

b'2 


U. f U 


0.40 


5.14 


6.24 


0.59 


0.35 


4.20 


K 1 A 

0.14 


S'3 


on 

0.80 


0.46 


6.23 


7.49 


0.68 


0.41 


5.06 


6.15 








af 


(po)(E^) (mb GeV 2 ) 








1 


2.60 


1.49 


19.01 


23.10 


2.23 


1.33 


15.64 


19.20 


Si 


2.25 


1.27 


16.03 


19.55 


1.94 


1.13 


13.30 


16.37 


s 2 


2.31 


1.24 


16.67 


20.22 


1.99 


1.11 


13.93 


17.03 


s 3 


2.61 


1.41 


20.22 


24.24 


2.24 


1.27 


16.78 


20.29 



TABLE IV. The minijet cross section, Eq. (]6|), first and second moments of the transverse 



energy distribution, Eqs. (17) and ( |19|) respectively with po = 2 GeV, integrated over b and r and 
divided by AB, within PHENIX, \y\ < 0.35. Note that the cross sections and moments are given 
over all azimuth. Results for both sets of parton distributions are separated into contributions from 
quarks, antiquarks and gluons as well as the total. The calculations are done without shadowing, 
5 = 1, and with shadowing parameterizations Si, S2, and S3. 
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GRV 94 LO 




MRST LO 








gluon 




total 


gluon 


total 




o 


no 


IHS 


HS 


irlo xlo 


IHS 


HS 












(6 = 0,po) (GeV/fm 3 ) 








1 


00c 


- 


920 


404 


- 


457 




c 

Si 


A AO 


413 


487 


454 219 


204 


248 


231 


c 

■>2 


48b 


457 


527 


496 246 


933 


973 


ZO 1 


S3 


566 


543 


620 


594 283 


272 


318 


306 








n{(b = 0,p ) (1/fm 3 ) 








1 


282 




309 


128 




144 




Si 


148 


138 


162 


151 68 


63 


77 


72 


So 


154 


145 


166 


155 73 


68 


81 


76 


S3 


183 


174 


200 


191 85 


82 


96 


92 



TABLE V. The energy density, Eq. (|2l|), and number density, Eq. (|22]), at b = from minijet 
production alone with po = 2 GeV within CMS, \y\ < 2.4 are given for both sets of parton 
distributions. Results are shown for homogeneous (HS) and inhomogeneous (IHS) shadowing, with 
the latter based on Sws- Both the gluon contribution alone and the total for gluons with three 
light quark flavors are presented. The calculations are done without shadowing, 5 = 1, and with 
shadowing parameterizations Si, S2, and S3. 



GRV 94 LO MRST LO 

gluon total gluon total 



s 


HS 


IHS 


HS 


IHS 


HS 


IHS 


HS 


IHS 










0,p ) (GeV/fm 3 ) 








1 


986 




1079 




426 




481 




Si 


513 


478 


561 


522 


226 


210 


256 


237 


So 


563 


531 


607 


572 


255 


240 


281 


264 


S3 


661 


634 


719 


689 


296 


285 


331 


317 








n((b 


= 0,p ) (1/fm 3 


) 








1 


332 




362 




134 




151 




Si 


170 


158 


186 


173 


70 


65 


79 


73 


s 2 


177 


167 


190 


178 


75 


70 


84 


77 


S3 


213 


203 


231 


220 


89 


85 


99 


95 



TABLE VI. The energy density, Eq. (^l]), and number density, Eq. (|22"|), at b = from minijet 
production alone with p$ = 2 GeV within ALICE, \y\ < 1 are given for both sets of parton 
distributions. Results are shown for homogeneous (HS) and inhomogeneous (IHS) shadowing, with 
the latter based on Sws- Both the gluon contribution alone and the total for gluons with three 
light quark flavors are presented. The calculations are done without shadowing, S = 1, and with 
shadowing parameterizations Si, S2, and S3. 
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GRV 94 LO 


MRST LO 








gluon 


total 


gluon 


total 






no 


IHS 


HS IHS HS 


IHS 


HS 










T 

e/(6 = 0,po) (GeV/fm 3 ) 








1 


ion 
18.9 


- 


23.2 - 15.7 


- 


19.5 




c 


1 k a. 
15.0 


15.3 


19.4 19.0 13.1 


12.8 


16.3 


16.0 


C 
■>2 


16.1 


1 5 9 


19 8 19^ 136 

j-iJ.(J .lo.u 


1 3 3 


1 6 8 


16.5 


S3 


19.6 


19.7 


23.8 23.9 16.5 


16.5 


20.2 


20.2 








n{(b = 0,p ) (1/fm 3 ) 








1 


7.2 




8.7 - 5.8 




7.2 




Si 


5.9 


5.8 


7.2 7.1 4.8 


4.7 


6.0 


5.8 


So 


6.0 


5.9 


7.3 7.2 5.0 


4.8 


6.1 


6.0 


S3 


7.3 


7.3 


8.8 8,8 6.0 


6.1 


7.4 


7.4 



TABLE VII. The energy density, Eq. (|21|), and number density, Eq. (22), at b = from minijet 
production alone with po = 2 GeV within the STAR calorimeter, |y| < 0.9 are given for both 
sets of parton distributions. Results are shown for homogeneous (HS) and inhomogeneous (IHS) 
shadowing, with the latter based on Sws- Both the gluon contribution alone and the total for 
gluons with three light quark flavors are presented. The calculations are done without shadowing, 
S = 1, and with shadowing parameterizations Si, S2, and S3. 







GRV 94 LO 




MRST LO 








gluon 


total 


gluon 


total 




s 


HS 


IHS 


HS 


IHS HS 


IHS 


HS 


IHS 










0,Po) (GeV/fm 3 ) 








1 


19.9 




24.3 


16.1 




19.9 




Si 


16.5 


16.2 


20.2 


19.8 13.5 


13.2 


16.6 


16.2 


So 


17.0 


16.8 


20.7 


20.3 14.0 


13.7 


17.1 


16.8 


S3 


20.7 


20.7 


24.9 


24.6 16.8 


16.9 


20.5 


20.6 








n((b 


= 0,p ) (1/fm 3 ) 








1 


7.6 




9.1 


6.0 




7.4 




Si 


6.2 


6.1 


7.5 


7.3 4.9 


4.8 


6.0 


5.9 


s 2 


6.3 


6.2 


7.6 


7.4 5.1 


5.0 


6.2 


6.1 


S3 


7.7 


7.7 


9.2 


9.2 6.2 


6.2 


7.5 


7.5 



TABLE VIII. The energy density, Eq. (^T|), and number density, Eq. (^2j), at b = from minijet 
production alone with pq = 2 GeV within PHENIX, |y| < 0.35 are given for both sets of parton 
distributions. Results are shown for homogeneous (HS) and inhomogeneous (IHS) shadowing, with 
the latter based on S\ys- Both the gluon contribution alone and the total for gluons with three 
light quark flavors are presented. The calculations are done without shadowing, 5 = 1, and with 
shadowing parameterizations S±, S2, and S3. 
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Detector Rapidity E? (mb GeV) 1| (mb GeV 2 ) 

CMS \y\ < 2.4 135 450 

ALICE \y\ < 1 56 188 

STAR \y\ < 0.9 34 112 

PHENIX \y\ < 0.35 13 44 



TABLE IX. The first and second Et moments of the soft contribution adjusted to the accep- 
tance of the experiments at the LHC and RHIC. We assume a p <? = 40 mb at RHIC and a^P = 60 
mb at the LHC. 



Detector a(S = 1) (nb) a(S = 5Q (nb) ajS = 5 2 ) (nb) ojS = 5 3 ) (nb) 

4 < m < 9 GeV 

CMS 4.05 1.90 1.57 2.26 

ALICE 1.89 0.86 0.68 1.04 

STAR 0.32 0.26 0.26 0.28 

PHENIX 0.13 0.10 0.10 0.11 

11< m < 20 GeV 

CMS 0.48 0.25 0.24 0.33 

ALICE 0.23 0.10 0.10 0.15 

TABLE X. Leading order Drell-Yan cross section, in units of nb per nucleon pair, integrated 
over all impact parameters, for the MRST LO parton densities. Full azimuthal coverage is assumed. 



Detector ajS = 1) (fib) a[S = 5Q (fib) a{S = S 2 ) (/ib) a{S = 5 3 ) (//b) 

Color Evaporation Model 

CMS 43.5 19.6 19.2 22.8 

ALICE 18.8 8.23 8.00 9.63 

STAR 1.62 1.12 1.12 1.43 

PHENIX 0.65 0.44 0.40 0.56 

NRQCD 

CMS 51.0 23.8 27.1 29.8 

ALICE 21.5 9.75 11.1 12.3 

STAR 1.54 1.11 1.16 1.49 

PHENIX 0.60 0.44 0.45 0.58 



TABLE XL J/tp production cross sections in the color evaporation and NRQCD approach in 
units of /ib per nucleon pair. No nuclear absorption of the J/tp in the final-state is included. 
The color evaporation cross sections were calculated with the MRST LO parton densities and the 
NRQCD results were obtained with the CTEQ 3L distributions. Both are normalized so as to 
agree with results from charmonium hadroproduction. Full azimuthal coverage is assumed. 
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Detector a{S = 1) (nb) a{S = 5Q (nb) ojS = S 2 ) (nb) a(S = 5 3 ) (nb) 

Color Evaporation Model 



CMS 


377 


187 


249 


267 


ALICE 


169 


80 


107 


117 


STAR 


4.80 


4.38 


4.72 


5.76 


PHENIX 


1.92 


1.77 


1.89 


2.36 






NRQCD 






CMS 


419 


282 


343 


365 


ALICE 


181 


119 


146 


157 


STAR 


6.19 


5.92 


6.17 


6.74 


PHENIX 


2.52 


2.43 


2.52 


2.78 



TABLE XII. T production cross sections in the color evaporation and NRQCD approach in 
units of /A) per nucleon pair. No nuclear absorption of the T in the final-state is included. The color 
evaporation cross sections were calculated with the MRST LO parton densities and the NRQCD 
results were obtained with the CTEQ 3L distributions. Both are normalized so as to agree with 
results from bottomonium hadroproduction. Full azimuthal coverage is assumed. 
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FIG. 2. The three shadowing parameterizations for A = 200 for (a) valence quarks, (b) sea 
quarks, and (c) gluons, relative to S = 1. The Si parameterization is shown in the solid curves. 
The S*2 ratios are given by the dashed curves. At low x, the lower curves are for Q = 2 GeV while 
the upper are for Q = 10 GeV. The S3 ratios, in the dot-dashed curves, are shown for the uy and 
u. The lower curves at low x are for Q = 1.5 GeV while the upper curves at low x are for Q = 10 
GeV. 
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FIG. 3. The radial distribution of shadowing for three models: SwS) Eq. @, (solid), Sr, Eq. (|5|) 
(dashed) and S p , Eq. (Q). All curves are normalized to a homogeneous S l (A, x, Q 2 ) of 0.7. Note 
that Sws is evaluated at z = 0. 
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FIG. 4. The rapidity distributions of quarks, antiquarks, gluons and the sum of all contributions 
in Pb+Pb collisions at ^/snn = 5.5 TeV integrated over b and divided by AB calculated with the 
GRV 94 LO parton distributions for po = 2 GeV. The solid curve is without shadowing, the dashed 
is with shadowing parameterization Si, the dot-dashed is with S*2 and the dotted uses S3. 
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FIG. 5. The rapidity distributions of quarks, antiquarks, gluons and the sum of all contributions 
in Pb+Pb collisions at ^/s^n = 5.5 TeV integrated over b and divided by AB calculated with the 
MRST LO parton distributions for po = 2 GeV. The solid curve is without shadowing, the dashed 
is with shadowing parameterization Si, the dot-dashed is with S*2 and the dotted uses S3. 
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FIG. 6. The rapidity distributions of quarks, antiquarks, gluons and the sum of all contributions 
in Au+Au collisions at ^/snn = 200 GeV integrated over b and divided by AB calculated with 
the GRV 94 LO parton distributions for p$ = 2 GeV. The solid curve is without shadowing, the 
dashed is with shadowing parameterization Si, the dot-dashed is with S*2 and the dotted uses S3. 
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FIG. 7. The rapidity distributions of quarks, antiquarks, gluons and the sum of all contributions 
in Au+Au collisions at ^Js^n = 200 GeV integrated over b and divided by AB calculated with the 
MRST LO parton distribution for po = 2 GeV. The solid curve is without shadowing, the dashed 
is with shadowing parameterization Si, the dot-dashed is with S2 and the dotted uses 6*3. 
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FIG. 8. The number of scatters suffered by an incoming gluon as a function of impact param- 
eter. The LHC results at x\ = KT 2 are shown in (a) and (b) for GRV 94 LO and MRST LO 
parton densities respectively while the RHIC calculations at X\ = 1CT 1 are shown in (c) and (d) 
for the GRV 94 LO and MRST LO sets. The solid curve is without shadowing, the dashed is with 
shadowing parameterization Si, the dot-dashed is with S2 and the dotted uses S3. 
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FIG. 9. The first E T moment, a(p )(E^), function of rapidity for quarks, antiquarks, 

gluons and the sum of all contributions in Pb+Pb collisions at ^/sjvjv = 5.5 TeV integrated over b 
and divided by AB calculated with the GRV 94 LO parton distributions for po = 2 GeV. The solid 
curve is without shadowing, the dashed is with shadowing parameterization Si, the dot-dashed is 
with S2 and the dotted uses S3. 
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FIG. 10. The first E T moment, a{p )(E J T ), function of rapidity for quarks, antiquarks, 

gluons and the sum of all contributions in Pb+Pb collisions at -^/sjvtv = 5.5 TeV integrated over b 
and divided by AB calculated with the MRST LO parton distributions for po = 2 GeV. The solid 
curve is without shadowing, the dashed is with shadowing parameterization Si, the dot-dashed is 
with S2 and the dotted uses S3 . 
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FIG. 11. The first E T moment, a{p )(E J T ), function of rapidity for quarks, antiquarks, 

gluons and the sum of all contributions in Au+Au collisions at ^snn = 200 GeV integrated over b 
and divided by AB calculated with the GRV 94 LO parton distributions for po = 2 GeV. The solid 
curve is without shadowing, the dashed is with shadowing parameterization Si, the dot-dashed is 
with S2 and the dotted uses S3. 
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FIG. 12. The first E T moment, a(p )(E J T ), function of rapidity for quarks, antiquarks, 

gluons and the sum of all contributions in Au+Au collisions at y/s^N = 200 GeV integrated over b 
and divided by AB calculated with the MRST LO parton distributions for po = 2 GeV. The solid 
curve is without shadowing, the dashed is with shadowing parameterization Si, the dot-dashed is 
with S2 and the dotted uses S3 . 
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FIG. 13. The impact parameter dependence of the first Et moment, a(po)(E^) relative to the 
E T moment with S = 1 in CMS, \y\ < 2.4, calculated with the MRST LO distributions with p = 2 
GeV. The upper plot shows the ratio for quarks and antiquarks, the middle plot is the gluon ratio 
and the lower plot is for the total. The horizontal lines show the homogeneous shadowing results: 
dashed for S\ , dot-dashed for S2 , and dotted for S3 . The inhomogeneous shadowing results for S\ , 
circles, S2, squares, and S3, diamonds are shown for S\ys (filled symbols) and S p (open symbols). 
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FIG. 14. The impact parameter dependence of the first Ex moment, a(jpo){E T ) relative to 
the E T moment with S = 1 in STAR, \y\ < 0.9, calculated with the MRST LO distributions with 
Po = 2 GeV. The upper plot shows the ratio for quarks and antiquarks, the middle plot is the gluon 
ratio and the lower plot is for the total. The horizontal lines show the homogeneous shadowing 
results: dashed for Si, dot-dashed for S2, and dotted for S3. The inhomogeneous shadowing results 
for Si, circles, S2, squares, and S3, diamonds are shown for Sws (filled symbols) and S p (open 
symbols) . 
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FIG. 15. The Et distribution predicted for the CMS detector in the interval \y\ < 2.4, calcu- 
lated with the MRST LO distributions and po = 2 GeV. The upper plot is for central collisions 
with b < 0.2Ra, the middle plot shows the region 0.9Ra < b < I.IRa, and the lower plot shows the 
entire Et distribution. The lines indicate the homogeneous shadowing results: solid for no shad- 
owing, dashed line for Si, dot-dashed for 52, and dotted for S3. The inhomogeneous shadowing 
results for Si, circles, 52, squares, and 53, diamonds are shown for 5ws- 
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FIG. 16. The .Et distribution predicted for the ALICE detector in the interval |y| < 1, calcu- 
lated with the MRST LO distributions and po = 2 GeV. The upper plot shows central collisions 
with b < 0.2Ra, the middle plot shows the region 0.9Ra < b < 1.1-R.4, and the lower plot shows the 
entire Et distribution. The lines indicate the homogeneous shadowing results: solid for no shad- 
owing, dashed line for Si, dot-dashed for 52, and dotted for S3. The inhomogeneous shadowing 
results for Si, circles, 52, squares, and 53, diamonds are shown for 5ws- 
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FIG. 17. The Et distribution predicted for the STAR detector in the interval \y\ < 0.9, calcu- 
lated with the MRST LO distributions and po = 2 GeV. The upper plot is for central collisions 
with b < 0.2Ra, the middle plot shows the region 0.9Ra < b < I.IRa, and the lower plot shows the 
entire Et distribution. The lines indicate the homogeneous shadowing results: solid for no shad- 
owing, dashed line for Si, dot-dashed for S2, and dotted for S3. The inhomogeneous shadowing 
results for Si, circles, 52, squares, and S3, diamonds are shown for Sws- 
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FIG. 18. The Et distribution predicted for the PHENIX detector in the interval \y\ < 0.35, 
calculated with the MRST LO distributions and po = 2 GeV. The upper plot emphasizes central 
collisions with b < 0.2Ra, the middle plot shows the region 0.9Ra < b < 1.1-Ra, and the lower 
plot shows the entire Et distribution. The lines indicate the homogeneous shadowing results: solid 
for no shadowing, dashed line for S±, dot-dashed for S2 , and dotted for S3. The inhomogeneous 
shadowing results for Si, circles, S2, squares, and S3, diamonds are shown for S\ys- 
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FIG. 22. The Drell-Yan rapidity distribution in the mass interval 4 < m < 9 GeV, relative to 
S = 1 for Pb+Pb collisions at the LHC, calculated with the MRST LO distributions. Central, 
b < 0.2Ra, semi-central, 0.9Ra < b < I.IRa, peripheral, 1.9Ra < b < 2.1Ra impact parameters 
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FIG. 23. The Drell-Yan rapidity distribution in the mass interval 4 < m < 9 GeV, compared 
to the distribution with S = 1 in Au+Au collisions at RHIC, calculated with the MRST LO dis- 
tributions. Central, b < 0.2Ra, semi-central, 0.9Ra < b < I.IRa, peripheral, 1.9Ra < b < 2.1Ra 
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neous shadowing result. The dashed line represents Si, the dot-dashed, S2, and the dotted, S3. 
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FIG. 35. The T rapidity distribution calculated in the NRQCD model with the CTEQ 3L 
distributions, compared to the distribution with S = 1 in Pb+Pb collisions at the LHC. Central, 
b < 0.2Ra, semi-central, 0.9Ra < b < I.IRa, peripheral, 1.9Ra < b < 2.1Ra impact parameters 
are shown along with the integral over all b. The lines show the homogeneous shadowing result: 
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